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We introduce the following family of evolutionary 1+1 PDEs that 
describe the balance between convection and stretching for small vis- 
cosity in the dynamics of ID nonlinear waves in fluids: 

nit + um x + bu x m = vm xx , with u = g * m . 



Here u = g*m denotes u(x) = g(x — y)m(y) dy . This convolution 
(or filtering) relates velocity u to momentum density m by integration 
against the kernel g{x). We shall choose g{x) to be an even function, 
so that u and m have the same parity under spatial reflection. When 
v = 0, this equation is both reversible in time and parity invariant. 
We shall study the effects of the balance parameter b and the kernel 
g{x) on the solitary wave structures, and investigate their interactions 
analytically for v = and numerically for small viscosity, v ^ 0. 

This family of equations admits the classic Burgers "ramps and 
cliffs" solutions which are stable for — 1 < b < 1 with small viscosity. 

For b < — 1, the Burgers ramps and cliffs are unstable. The stable 
solution for b < — 1 moves leftward instead of rightward and tends to 
a stationary profile. When m = u — a 2 u xx and v = 0, this profile is 



email: dholm@lanl.gov 



February 26, 2002 



Abstract 




convection stretching viscosity 



1 



given by u(x) ~ sech 2 (x/(2a)) for b = —2, and by u(x) ~ sech(x/a) 
for b = —3. 

For b > 1, the Burgers ramps and cliffs are again unstable. The 
stable solitary traveling wave for b > 1 and v = is the "pulson" 
u(x, t) = cg(x—ct), which restricts to the "peakon" solution in the spe- 
cial case g(x) = e~\ x \l a when m = u — a 2 u xx . Nonlinear interactions 
among these pulsons or peakons are governed by the superposition of 
solutions for b > 1 and u = 0, 

N 

u(x,t) = ^2pi(t)g(x - qi(t)) . 
i=i 

These solutions obey a finite dimensional dynamical system for the 
time-dependent speeds pi(t) and positions qi(t). We study the pul- 
son and peakon interactions analytically, and we determine their fate 
numerically under adding viscosity. 
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1 Introduction 

1.1 The b-family of fluid transport equations 

We shall analyze a one- dimensional version of active fluid transport that is 
described by the following family of 1+1 evolutionary equations, 

m t + um x + bu x m = 0, with u = g*m, (1) 
convection stretching 

in independent variables time t and one spatial coordinate x. 

We shall seek solutions for the fluid velocity u(x,t) that are defined ei- 
ther on the real line and vanishing at spatial infinity, or on a periodic one- 
dimensional domain. Here u = g *m denotes the convolution (or filtering), 

/oo 
g(x -y)m(y)dy, (2) 
-oo 

which relates velocity u to momentum density m by integration against kernel 
g(x) over the real line. We shall choose g(x) to be an even function, so that 
u and m have the same parity. 

The family of equations ([!]) is characterized by the kernel g and the real 
dimensionless constant b, which is the ratio of stretching to convective trans- 
port. As we shall see, b is also the number of covariant dimensions associated 
with the momentum density m. The function g(x) will determine the trav- 
eling wave shape and length scale for equation ([]]), while the constant b will 
provide a balance or bifurcation parameter for the nonlinear solution be- 
havior. Special values of b will include the first few positive and negative 
integers. 

The quadratic terms in equation ([!]) represent the competition, or bal- 
ance, in fluid convection between nonlinear transport and amplification due 
to b— dimensional stretching. For example, if m is fluid momentum (a one- 
form density in one dimension) then b = 2. Equation ([!]) with b = 2 arises 
in the nonlinear dynamics of shallow water waves, as shown in j| and @. 
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Equation (j^) with 6 = 2 and 6 = 3 appears in the theory of integrable partial 
differential equations 0, || ||. The three-dimensional analog of equation ([[J) 
with 6 = 2 was introduced in || [J. Applying the proper viscosity to this 
three-dimensional analog with 6 = 2 produces the Navier-Stokes-alpha model 
of turbulence f3j. The ID version of this turbulence model is 

m t + um x + bu x m = vm xx , with u = g * m . (3) 
convection stretching viscosity 

We shall compare our analysis of equation (JJ) with numerical simulations of 
for small viscosity. 

1.2 Outline of the paper 

After summarizing previous investigations of particular cases in the b-family 
(jl|) of active transport equations, section || discusses its symmetries and other 
general properties such as parity and reversibility. Section discusses the 
traveling waves of equation (|l[) and derives their Pulson solutions, which 
may be generalized functions for 6 > 1. Section [| analyzes the interaction 
dynamics of the Pulson solutions for any positive 6 and any g. Section 
[5] specializes the analysis of the Pulson solutions to the Peakons, for which 
g(x) = e~\ x \l a is a peaked pulse of width a, and 6 is taken to be arbitrary. In 
section ^| we add viscosity to the peakon equation, and describe our numerical 
methods for illustrating the different types of behavior that may arise in the 
initial value problems for Peakon solutions with 6 > 0, 6 = and 6 < — 1. 
Section using these numerical methods to determine how viscosity affects 
the fate of the peakons. Section § provides a synopsis of the figures. Section 
[| summarizes the paper's main conclusions. 

2 History and general properties of the b- 
equation 

Camassa and Holm derived the following equation for unidirectional mo- 
tion of shallow water waves in a particular Galilean frame, 

m t + um x + 2 u x m = — cqu x — 7 u xxx & m = u — a 2 u xx . (4) 
convection stretching dispersion 
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Here m = u — a 2 u xx is a momentum variable, partial derivatives are denoted 
by subscripts, the constants a 2 and 7/c are squares of length scales, and 
c o = VFh is the linear wave speed for undisturbed water of depth h at rest 
under gravity g' at spatial infinity, where u and m are taken to vanish. Any 
constant value u = u is also a solution of (U). 

Equation (^) was derived using Hamiltonian methods in |2[ and was shown 
in || also to appear as a water wave equation at quadratic order in the 
standard asymptotic expansion for shallow water waves in terms of their two 
small parameters (aspect ratio and wave height). The famous Korteweg-de 
Vries (KdV) equation appears at linear order in this asymptotic expansion 
and is recovered from equation ([|) when a 2 — > 0. Both KdV at linear 
order and its nonlocal, nonlinear generalization in equation (^) at quadratic 
order in this expansion have the remarkable property of being completely 
integrable by the isospectral transform (1ST) method. The 1ST properties 
of KdV solitons are well known and these properties for equation (|j) were 
studied in, e.g., and Q. 

When linear dispersion is absorbed by a Galilean transformation and 
a velocity shift, equation @ reduces to an active transport equation that 
contains competing quadratically nonlinear terms representing convection 
and stretching, 

m t + um x + 2 u x m = , with m = u — a 2 u xx . (5) 
convection stretching 

This is a special case of equation ([I]) for which 6 = 2 and g(x) = e~\ x \l a . The 
traveling wave solution of @ is the "peakon," u(x,t) = ce -\ x ~ ct \l a found 
in 0, where e~\ x \l a is the Green's function for the Helmholtz operator that 
relates m and u. The interactions among iV peakons are governed by the 
2iV dimensional dynamical system for the speeds Pi(t) and positions q^t) 
i = 1, . . . , N, appearing in the solution, 

N 

u(x,t) = J2Pi(t)e~ lx ~ 9im - (6) 

i=l 

As shown in Camassa and Holm , a closed integrable Hamiltonian system of 
ordinary differential equations for the speeds Pi(t) and positions gj(i) results 
upon substituting the superposition of peakons (0) into equation (|3p. This 
integrable system governs the dynamics of the peakon interactions. 
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A variant of equation (|^) with coefficient 6 = 2 — ^6 = 3, 

m t + um x + 3u x m = 0, with m = u — a 2 u xx , (7) 
convection stretching 

was first singled-out for further analysis by Degasperis and Procesi ||. De- 
gasperis, Holm and Hone Q discovered that this 2 — > 3 variant of equation (g) 
possesses superposed peakon solutions (|6]) and also is completely integrable 
by the isospectral transform method. Thus, the N— dimensional pulson solu- 
tion (H) is a completely integrable dynamical system for (0), as well, but with 
different dynamics for the speeds pi(t) and positions qi(t) of the peakons. 

Fringer and Holm extended the zero-dispersion equation (|5]) for the 
peakons to the "pulson" equation, 

m t + um x + 2 u x m = , with u = g * m . (8) 
convection stretching 

Here u = g * m denotes the convolution (or filtering) 

/■oo 

u(x) = / g(x - y)m(y) dy , (9) 



that relates velocity u to momentum density m by integration against the 
kernel g(x). Fringer and Holm chose g(x) to be an even function, so that 
u and m have the same parity. They studied the effects of the shape of the 
traveling wave u(x,t) = cg(x — ct) on its interactions with other traveling 
waves in the superposed solution, 

N 

u{x,t) = ^vti) g{x - <n{t)) , (io) 

i=l 

This superposed solution of traveling wave forms with time dependent speeds 
Pi(t) and positions qi(t), i = 1, . . . , N, revealed that the nonlinear interactions 
among these pulsons occur by elastic two-pulson scattering, even though the 
Fringer-Holm pulson equation (H) is not integrable for an arbitrary choice of 
the kernel g. When g(x) = e~^ x ^ a is assumed, the pulson equation for 6 = 2 
in (H) specializes to the peakon equation for b = 2 in (|). 
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2.1 Discrete symmetries: reversibility, parity and sig- 
nature 

Equation ([]]) for m is reversible, or invariant under t — *■ — t and u — > —u. 
The latter implies m — > —m. Hence, the transformation u(x,t) — > —u(x, —t) 
takes solutions into solutions, and in particular, it reverses the direction and 
amplitude of the traveling wave u(x, t) = cg(x — ct). 

We choose g (x) to be an even function so that m and u = g*m both have 
odd parity under mirror reflections. Hence, equation ([I]) is invariant under 
the parity reflections u(x, t) — > — u(—x, t), and the solutions of even and odd 
parity form invariant subspaces. 

Equation (|l]) implies a similar reversible, parity invariant equation for the 
absolute value \m\, 

\m\t + u\m\ x + bu x \m\ = , with u = g * m . (11) 

So, the positive and negative components m± = |(m ± |m|) satisfy equation 
([]]) separately. Also, if m is initially zero, it remains so. This is conservation 
of the signature of m. 

2.2 Lagrangian representation 

If m 1//fe is well-defined, equation (|l]) may be written as the conservation law 

dtm 1 ^ + d x {m^ b u) = 0, (12) 
and equation ([TT|) for the absolute value implies 

d t \m\ 1/b + d x {\m\ 1/b u) = 0. (13) 
Adding and subtracting equations ( |T2"D and flllf ) implies 

d t (m 1/6 )± + d x ((m 1/b )±u) = with {m l ' h ) ± = ^(m 1 / 6 ± \m\^ b ) . (14) 

Consequently, regions of positive and negative m are transported by the same 
velocity and their boundaries propagate so as to separately preserve the two 
integrals, 

/oo 
(m 1 ' b ) ± dx. (15) 
-oo 
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The common transport velocity allows a transformation to Lagrangian coor- 
dinates X± defined by 



dX± = (m 1/b )±(dx - udt) so that d t X± + ud x X± = . (16) 

This formal transformation is not strictly defined where (m 1 / b )± vanishes. 
However, by equation (0), regions where (m 1 ^ b )± vanishes do not propagate 
and do not contribute to the integrated value of X± = J^ 00 (m 1 ^ b )±(y, 0) dy. 
Hence, these regions may be identified and excluded initially. The formal 
inverse relation holding in the remaining regions, 

d x = (m 1/b )± x dX± + udt , (17) 

implies that 

rlnr 

u(x,t), (18) 



dx 
~dt 



so the Lagrangian trajectories x = x(X±,t) of positive and negative inte- 
grated initial values of X± = f^ 00 (rn 1 ^ b )±(y, 0) dy are transported by the 
same velocity u = g * m. 

2.3 Preservation of the norm ||m||jri/& for < b < 1 

If |m| 1//fc is well-defined, the continuity equation form (|]|) of equation ([!]) 
implies conservation of 

l 1 ^ dx = / | mo | dx , where m (x) = m(x, 0) . (19) 



oo poo 

\m\ 

■oo J — oo 



This integral is conserved for all 6, but only defines a norm (the L x l b norm 
Umll^i/b) in the closed interval < b < 1. In the limit b — > this becomes 
the norm |m| max . Hence, when 6 = equation ([!]) has both a maximum 
principle and a minimum principle for m. Such a principle is meaningful only 
if m 1//fc is an ordinary function, e.g., if m is not a generalized function, such 
as the delta functions that occur for the peakons we will discuss below. 

Thus, the L l / b norm is conserved by equation (|]) provided Imp b 

is well-defined for the closed interval < b < 1. One may also define the 
corresponding conserved norm for 1/m in the closed interval — 1 < b < 0, 
provided Im] 1 / 6 is well-defined on this interval. 
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2.4 Lagrangian representation for integer b 

Fluid convection means transport of a quantity by the fluid motion. Exam- 
ples of transported fluid quantities are circulation (a one-form) in Kelvin's 
theorem for the Euler equations and its exterior derivative the vorticity (a 
two-form, by Stokes theorem) in the Helmholtz equation. For a Lagrangian 
trajectory x(X,t) satisfying x(X, 0) = X and 

dx = {m l/b )- l dX + udt , (20) 

we have seen that the conservation law (|T2"D implies 

m 1/b (x,t)dx = m l/b (X,0)dX, (21) 

provided that m 1 ^ is a well defined function. The last issue may be circum- 
scribed as follows when b is an integer. In ID, higher order differential forms 
may be created by using the direct, or tensor, product, e.g., dx®dx = dx® 2 . 
Consequently, the tensor product of each side of equation (^1|) b times givesQ 

m(x, t)dx® b = m(X, 0)dX® b . (22) 

Taking the partial time derivative of this equation at constant Lagrangian 
coordinate X and using dx/dt\ x = u yields equation ([I]) in the form 

(m(x, t)dx® b ^J = (m t + um x + b u x m)dx® b = . (23) 



dt 



x 



Thus, when the parameter b in equation (Q) is an integer, it may be regarded 
geometrically as the number of dimensions that are brought into play by 
coordinate transformations of the quantity mdx® b associated with m. Cases 
of equation (0) with negative integer b < may be interpreted as 

d / m(d x r { - b) )=0. (24) 



(md x ) = (m t + um x - u x m) (d x ) = , (25) 



dt x 

For example, the case b = — 1 may be written as 
d 
dl 

in which the difference of terms (um x — u x m)d x is the commutator of the 
vector fields ud x and md x on the real line. The rest of the paper will remain 
in the Eulerian (spatial) representation. 

1 Cases with integer values of b will allow to to be a generalized function. Cases with 



non-integer values of b will revert to equation (21) for which to is required to be a classical 
function. 



10 



2.5 Reversibility and Galilean covariance 

Equation ([D is reversible, i.e., it is invariant under the discrete transforma- 
tion u(x,t) — > —u(x,—t). Equation ([!]) is also Galilean-covariant for all 6. 
In fact, equation ([[]) keeps its form under transformations to an arbitrarily 
moving reference frame for all b. This includes covariance under transforming 
to a uniformly moving Galilean frame. However, only in the case b = is 
equation (H) Galilean invariant, assuming that m Galileo-transforms in the 
same way as u. In this case, equation (JJ) transforms under 

t — > t + to , x — > x + xq + ct , u — > u + c + Mo , m — > m + c + uo , (26) 

to the form 

nit + um x + bu x m + Mo^z + bu x (c + «o) = , with u = g * m . (27) 

Thus, equation (P is invariant under space and time translations (constants 
Xq and to); covariant under Galilean transforms (constant c), and acquires 
linear dispersion terms under velocity shifts (constant u ). Equation (|l]) 
regains Galilean invariance if m is Galilean invariant. However, the disper- 
sive term u m x introduced by the constant velocity shift u ^ breaks the 
reversibility of equation (JJ) even if m is invariant under this shift. 



2.6 Integral momentum conservation 

Equation ([12]) implies that M = m dx is conserved for any g when 6 = 1. 
However, when g(x) is even, the family of equations ([TJ) also conserves the 
total momentum integral M for any b. This is shown by directly calculating 
from (HD that 

j /»oo /»oo /*oo 

— / m(x) drr = (1 — b) / / m(x) g ' (x — y)m(y) dx dy = , (28) 



dl _ 

in which the double integral vanishes as the product of an even function and 
an odd function under interchange of x and y, when g '(—x) = —g '(x) . Hence, 
for even g(x), M = mdx is conserved for either periodic or vanishing 
boundary conditions and for any b. We shall assume henceforth that g(x) is 
even and, moreover, that the integral J m g * mdx is sign-definite, so that it 
defines a norm, 

/•oo roo /*oo 

||m||g = / mg*mdx= / m(x) g(x — y) m(y) dx dy > . (29) 



oo J — oo 



This norm is conserved by equation ([!]) when 6 = 2. 
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3 Traveling waves and generalized functions 

Its invariance under space and time translations ensures that equation ([I]) 
admits traveling wave solutions for any 6. Let us write the traveling wave 
solutions as 

u = u(z) and m = m(z) , where z = x — ct , (30) 
and let prime ' denote d/dz. 

3.1 Case 6 = 

3.1.1 Pulsons for 6 = 

For 6 = 0, equation (P is Galilean invariant and its traveling wave solutions 
satisfy 

(u(z) — c)m'(z) = , z = x — ct , (31) 



where prime ' denotes d/dz. Equation (|3TD admits generalized functions 
m'(z) ~ 5(z) matched by u — c = at z = 0. The velocity u is given by the 
integral of the Green's function that relates m and u = g * m, 



u — c ~ c 



g(y)dy . (32) 



3.1.2 Peakons, ramps, and cliffs for 6 = 

When g(x) = e^ x ^ a (the Green's function for the ID Helmholtz operator), 
we have m = u — a 2 u xx . Consequently, the equation m' = u' — a 2 u"' = ±26 (z) 
with u — c = 0at^ = 0is satisfied by 



u — c = ±c 



J e-\ y \ /a dyj * = ±c sgn(z) (l - e - |z|/a ) . (33) 



This represents a rightward moving traveling wave that connects the left 
states u — c = ±c to the same two right states. 

Definition 3.1 (Peakons) The symmetric connections u = ±ce~' z '/ a ; with 
a jump in derivative at z = 0, are the peakons, for which m = u — a 2 u xx and 
g(x) = e~' x ''/ a . 
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Definition 3.2 (Cliff's) The antisymmetric connections u = ±csgn(z)(l — 
e -l z l/ Q ) (withu — c = ±c connecting tou — c = ^c), with no jump in derivative 
at z = 0, are the regularized shocks (cliff's). These propagate rightward but 
may face either leftward or rightward, because equation (j^) in the absence 
of viscosity has no entropy condition that would distinguish between leftward 
and rightward facing solutions. 

Definition 3.3 (Ramps) Equation also has ramp-like similarity solu- 
tions u ~ x/t when g(x) = e~^ x ^ a for any b. These may emerge in the 
initial value problem for the peakon case of equation ([Jj and interact with 
the peakons and cliffs. 

Remark 3.4 (First integral for b = traveling waves) For 6 = 0, the 

traveling wave equation ( jffil j apparently has only the first integral for m = 

2 

Ti Gt Uj£2£ , 

(u -c)(u- a 2 u") - — + —v! 2 = K . (34) 

Thus, perhaps surprisingly, we have been unable to find a second integral for 
the traveling wave equation for peakons when b = 0. 

Remark 3.5 (Reversibility) Reversibility means that equation (|7p is in- 
variant under the transformation u(x,t) — > —u(x,—t). Consequently, the 
rightward traveling waves have leftward moving counterparts under the sym- 
metry c — > — c. The case of constant velocity u = ±c is also a solution. 

Figure [I] shows that the ramp and cliff pattern develops in the veloc- 
ity profile under the peakon equation ([I]) with g(x) = e~' x '/ a for a set of 
Gaussian initial conditions (Sy^) -1 exp(— (x — 50) 2 /w) of increasing widths 
w = 2.5, 5, 10, for a = 1 and 6 = 0. Apparently, the ramp solution is numer- 
ically stable, but the coexisting peakon solution is not stable in this case. A 
complete stability analysis of these various solutions is outside the scope of 
the present paper. Instead we shall investigate the solutions of equation ([[J) 
by numerically integrating selected examples. 
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Figure 1: Ramps and cliffs for 6 = 0. Inviscid 6-family, b = 0, a = 1, 
initial width w = 2.5, 5, 10. 



3.2 Case 6^0 

For b 0, the conservation law ( |12|) for traveling waves becomes, 



((u - c)m}l b ^ ' = , (35) 

which yields after one integration 

(u-c) b m = K, (36) 

where K is the first integral. For g(x) = e~^ x ^ a , so that m = u — a 2 u xx , this 
becomes 

{u-c)\u-a 2 u")=K . (37) 
For n - c ^ we rewrite this as 

a 2 u" = u- K(u-c)- b (38) 
and integrate again to give the second integral in two separate cases, 



a u 



2„./2_ I u 2 -^ E (u-c) 1 - b + 2H for 6^1, 
u 2 -2K\og(u- c) + 2H for 6 = 1. 
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We shall rearrange this into quadratures: 
dz du 



± 



a 



u 2 - ^- b ( u - C y- b + 2H 



1/2 



for b^l 



(39) 



and 



± 



dz 



a 



du 

■ 71/2 
u 2 - 2K\og(u - c) + 2H 



for 6=1 



(40) 



For 6 = 1 and K 7^ 0, the integral in equation (|4"0| ) is transcendental. 

3.2.1 Special cases of traveling waves for b ^ 

• For K = the two quadratures (|39D and ( [10|) are equal, independent 
of b, and elementary, thereby yielding the traveling wave solutions 



\z\/a 



U 



+ Vu 2 + 2H 



c + Vc 2 + 2H 



(41) 



with u — c = 0&tz = 0. 



For H = equation (O) recovers the peakon traveling wave. 



For H > equation (41) gives a rightward moving traveling wave that 
is a continuous deformation of the peakon. 

For H > and c = equation ( |4"1"D gives stationary solutions of the 
form 



+ vV + 2H ~ e" |z|/a . 



(42) 



3.3 Case 6 > 
3.3.1 Pulsons for 6 > 

Equation (|I]) for 6 > has nontrivial solutions vanishing as \z\ 
allow K = in equation (|3^), so that 



00 that 



(w — c) b m = . 



(43) 
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This admits the generalized function solutions 



m = c5(z) and u = g * m = cg(z) , (44) 

matched by u — c = at z = 0. This is the pulson traveling wave, whose 
shape in u is given by the kernel g. The constant velocity case u = c is a 
trivial traveling wave. 

Remark 3.6 (Pulson and peakon traveling waves) The pulson solution 



requires g(0) = 1 and g'(0) = 0. We shall assume for definiteness that 
the even function g(z) achieves its maximum at g(0) = 1, so that the sym- 
metric pulson traveling wave u(x,t) = cg(x — ct) moves at the speed of its 
maximum, which occurs at its center of symmetry. For example, the peakon 
u(x,t) = ce~\ x ~ ct \ moves at the speed of its peak. 

3.3.2 Peakons for b > 1 

Equation fl37|), for which g(x) = e~^ x '' a , yields the peakon traveling wave 

u(z) = ce" |2|/a and m(z) =u-a 2 u" = 2cS(z/a) (45) 
when K = 0. 

Figures and |3] show the effects of varying the width to of a Gaussian 
initial condition (Sy^) -1 exp(— (x — 50) 2 /w) for the peakon equation in a 
periodic domain, when a = 1 and b = 2, 3. As the width of the initial 
Gaussian increases, the figures show that more peakons of width a = 1 are 
emitted. (This is consistent with conservation of momentum.) The peakons 
are observed to be stable for b > 1, they propagate as solitary traveling 
waves, and they interact elastically. We shall discuss the peakon interactions 
in more detail in sections | through 

3.4 Case b < 

We shall examine the cases b = —0.5, —1, —2, —3, —4. Numerical results for 
b = —2 and b = —3 are described in section |3.4.6 . For other values of 



b < the analysis is similar, but it involves less elementary considerations 
such as transcendental or hyperelliptic functions. The numerics shown later 
will demonstrate that the elementary solutions discussed here, many of them 
stationary, do tend to emerge in numerical integrations of the initial value 
problem for equation (P with b < — 1. 
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u(x), w=10 




Figure 2: Peakons for 6 = 2. Inviscid 6-family, b — 2, a — 1, initial width 
w = 2.5,5,10. 
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u(x), w-2.5 



u(x), w=5 



u(x), w=10 




Figure 3: Peakons for 6 = 3. Inviscid 6-family, b — 3, a — 1, initial width 
w = 2.5,5,10. 
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U(X) b=-1/2 W=1 0,1 5,20 1=0.500.1000,1500,2000.2500 

0.15, ■ 1 i -i 




X 

Figure 4: Ramps and cliffs for b = —1/2. Inviscid 6-family, b = —1/2, 
a = 1, initial width w = 10, 15, 20. 



3.4.1 Case b = -1/2 

Figure [| shows that a ramp and cliff pattern develops in the velocity profile 
under the peakon equation ([I]) with g(x) = e~\' x ^ a for a set of Gaussian initial 
conditions (Sa/vt)" 1 exp(— (x — 100) 2 /w) of increasing width w = 10, 15,20, 
for a = 1 and b = —1/2. Apparently, the ramp solution is numerically stable 
for b = -1/2. 
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3.4.2 Case b = -1 

For b = — 1, equation fl39|) becomes 



^ = 7—7 — TIT.- < 46 > 



w 2 - K{u - c) 2 + 2H 



which integrates to 



-|*|/a _ « 



+ v/m 2 - ^ - c) 2 + 2H + Kc 



c + Vc 2 + 2# + 



(47) 



with it — c = at 2 = 0. (K = and H = recovers the peakon traveling 
wave.) 

Remark 3.7 (Stationary plane wave solutions for b = — 1) Equation (Q) 
/or 6 = — 1 is satisfied for any wavenumber k by, 

m = cos(k(x — ct) + 4>o) and u = g{k) cos{k(x — ct) + <f> ) , (48) 

where g(k) is the Fourier transform of the kernel g(x) and <po is a constant 
phase shift. In the absence of linear dispersion, these solutions are stationary, 
c = 0. When linear dispersion is added to equation ([][), these solutions are 
the ID analogs of Rossby waves in the 2D quasigeostrophic equations. 

Figure |5| shows the velocity profiles under evolution by the peakon equa- 
tion, ([!]) with g(x) = e~\ x " a , for a set of Gaussian initial conditions of in- 
creasing width w = 10, 15, 20 for a = 1 and b = —1. Evidently, the coexisting 
peakon solution for b = — 1 does not emerge because K ^ and H ^ for 
this initial condition. Instead, the stable solution is essentially stationary 
with a slight rightward drift and leaning slightly to the right. The reason for 
this lethargic propagation becomes clear upon writing the b-equation solely 
in terms of the velocity u(x, t) as 

u t + {b + l)uu x = a 2 (u xxt + uu xxx + bu x u xx ) (49) 

,„ / b — 1 2 
= a d x [u xt + uu xx H — u x 

6-3 



2 o2 / ^ 2 

a o x \u t + uu x H — u x 
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Figure 5: Stationary solutions for 6 = — 1. Inviscid 6-family, 6 = — 1, 
a = 1, initial width 10 = 10, 15, 20. 



Remark 3.8 (6 = —1 is a turning point) When b = — 1 £/ie nonlinear 
steepening term {b + l)uu x vanishes in and the residual propagation 
is due to its nonlinear "curvature terms" with higher order derivatives. In 
the parameter regime b > — 1 (resp. b < —1) the solutions of equation (W 
or ffiQ ) move rightward (resp. leftward), provided the curvature terms on the 
right hand side of equation ( $ty) are either negative, or sufficiently small. 



Remark 3.9 (Short wave limit equation) The high wavenumber, or short 
wave, limit of equation is 



dl[u t + uu x 



b-3 



0. 



(50) 



For 6 = 2 and linxj^-oo u x = 0, this integrates to become the Hunter-Saxton 
equation [WJ . For 6 = 3, it is the second derivative of the Burgers equation. 



3.4.3 Case 6 = — 2 stationary solutions 

For 6 = —2, the traveling wave quadrature ( p9|) becomes an elliptic integral 

dz d{u — c) 



± 



1/2 • 



(51) 
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u(x). w=10 



u(x), w=15 
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Figure 6: Stationary solutions for b = — 2. Inviscid 6-family, b = —2, 
a — 1, initial width w = 10, 15, 20. 



The hyperbolic limit of this equation for H = vanishes at infinity for the 
stationary solution (c = 0) to give 



u(z) = -^-sech 2 — 
K J 2K 2a 



(52) 



3.4.4 Case b = — 3 stationary solutions 

For b = —3, the hyperbolic limit H = of equation 

dz du 



is 



± 



a 



u 2 -f {u-cY 



1/2 



(53) 
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which for c = is 



± 



dz 



a 



u 



du 



1/2 ' 



(54) 



and may be integrated in closed form to obtain a continuous deformation of 
the peakon, 



-MA* 



u. 



1 + K/2 1 + y/l- Ku 2 /2 

Rearranging equation (|5^) and scaling u by uq gives 

Wo 



for b = -3 , c = and H = 0. (55) 



u(z) 



A e \z\/a + (! _ A) e -\z\/a 



(56) 



with A = 1 + y/\ - A/2, so that A G [1, 2] for A" e [0, 2}. For A = 1, we 
have u(z) = u sech (z/at). And for A = 2, we recover the stationary peakon, 



■urn = uoe 



\z\/a 



3.4.5 Case 6 = —4 stationary solutions 

For 6 < —4 the analytical expressions for the cnoidal waves become less ele- 
mentary, because the integral in equation (^) is then hyperelliptic. However, 
our numerics show that the dynamical behavior for b = —4 is similar to that 
of the cases b = —2 and b = —3 shown in Figures ^-0. Namely, a series of 
transient leftward propagating pulses, or leftons, of width alpha emerge and 
tend to a nearly steady state. Consistent with momentum (area) conserva- 
tion and the tendency toward pulses of width alpha, the number of emerging 
leftons increases with the width of the initial Gaussian. At a longer time 
scale, this train of pulses appears to tend toward stationary (c — >• 0). 



3.4.6 Numerical Results for b = —2 and b = — 3 

Figures ^ and |7] show that a series of leftons in the velocity profile emerges 
under the peakon equation for a set of Gaussian initial conditions of increas- 
ing width w = 10, 15, 20, for a = 1 and b = —2, —3. Apparently these are 
not peakons, because the velocity at which they move is not equal to their 
height. The leftons emerge from the initial Gaussian in order of height and 
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u(x). w=10 



u(x), w=15 



u(x), w=20 




Figure 7: Stationary solutions for b = — 3. Inviscid 6-family, b = — 3, 
a — 1, initial width w = 10, 15, 20. 
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then tend toward a nearly stationary state. The number of emerging pulses 
increases with the width of the initial Gaussian, as expected from momentum 
(area) conservation and the tendency toward pulses of width alpha, and the 
leftward speed of the emerging pulses increases with the magnitude of b. The 
latter is consistent with the coefficient (6+ 1) of the nonlinearity in equation 
( |T9| ) as b becomes more negative. 

Figure || shows the leftons at time T = 2500, versus u(x) ~ sech 2 (x/(2a)) 
for b = —2, and versus u(x) ~ sech(a;/a) for b = —3. By this time, the leftons 
have become stationary solutions with c = for both b = —2 and b = —3. 



4 Pulson interactions for b > 

As we have seen in section 3.3. 1| , the b-family of equations (P admits the 
pulson traveling wave solution ([PJ) for b > 0. The interaction dynamics 
among N of these pulsons is obtained by superposing the traveling wave 
solutions u(x, t) = cg{x — ct) as 

N N 

u(x, t) = y^ j p i {t)g{x - qj{t)) and m(x,t) = ^2pi(t)5(x - q^t)) , (57) 
i=i i=i 

for any b > and u = g * m, where the function g is even so that g'(0) = 
and bounded and we may set g(0) = 1. For these superpositions of pulsons 
to be exact solutions, the time dependent parameters Pi(t) and qi(t) must 
satisfy the following iV— dimensional particle dynamics equations obtained 
by substituting ( p7[) into equation ([[]), 



Pi 



N ar 

(1 - b) p^Pj 9 '(« - Hi) = (1 - b) , (58) 
i=i Qi 

YjPj9{^-qj) = ^-- (59) 



Here the generating function Gn is obtained by restricting the norm ||m|| 2 
in (p9|) to the class of superposed traveling wave solutions (0), as 

1 N 

Gn = - PiP i 9{Qi ~Qj)- (6°) 
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For b=-2: computed u(x) vs. u(x)~ sech 2 (x/(2a)) 
0.4i 1 " 1 1 1 
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For b=-3: computed u(x) vs. u(x)~ sech(x/a) 
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Figure 8: Stationary solutions for b = — 2 versus u(x) ~ sech 2 (a;/(2a)), 
and for b = — 3 versus u(x) ~ sech(x/o;). Inviscid Warmly, b = —2,-3, 
a = 1, initial width w = 10. 
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Thus, the symmetric kernel g(x) determines the shape of the traveling wave 
solutions (57D, an d these traveling waves interact nonlinearly via the pulson 
dynamics of pi(t) and with i — 1, . . . , N in equations fl58|) and fl59|) for 
b > 0. We shall see that the character of these interactions depends vitally 
on the value of b. 



4.1 Pulson interactions for 6 = 2 

When 5 = 2, equations and (^) describe the canonical dynamics of 
a Hamiltonian system with N degrees of freedom. These are the geodesic 
pulson equations studied in Fringer and Holm 0, in which the following 
results are obtained, 

• Equation ([TJ) conserves the kinetic energy = | J^rng * mdx. 

• Equations (^) and (^) describe canonical geodesic motion in an N— dimensional 

configuration space whose co- metric is g^(q) = g(qi — qj). 

• The generating function Gn is the kinetic energy Hamiltonian for the 

canonical geodesic motion. 

• The solutions in ([)?]) behave as particle-like pulses whose pairwise inter- 

actions as determined by equations (^) and ( [5§D comprise nonlinear, 
but elastic, scattering events. 

• The pairwise interactions for the pulsons can be solved analytically for 

any symmetric function g(x). 

Remark 4.1 As we shall show, the last two statements also hold for any 
b>l. 



4.2 Peakon interactions for 6 = 2 and 6 = 3: numerical 
results 

• Figure |9] shows the evolution of the velocity profiles in the 2-peakon 
and 3-peakon interactions for b = 2, with g(x) = e~^ x ^ a and a pe- 
riodic domain. The 3-peakon interaction decomposes into a series of 
2-peakon interactions. These simulations verify the analytical results 
for the 2-peakon interaction to three significant figures over propagation 
distances of about sixty peakon widths. 
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Figure 9: Peakons of width a for b = 2: collisions. Inviscid 6-family, 
b = 2, a = 5, initial width w — 5. 
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Figure 10: Peakons of width a for 6 = 3: collisions. Inviscid 6-family, 
b = 3, a = 5, initial width w — 5. 
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u(x), b=2 



u(x), b=3 




Figure 11: Peakons of width > a for b = 2,3: emergence of width-a 

peakons. Inviscid 6-family, b = 2, 3, a = 1, initial width w — 5. 



• Figure shows the evolution of the velocity profiles in the 2-peakon 
and 3-peakon interactions for 6 = 3, with g(x) = e~' x '/ a and a periodic 
domain. 

• Figure |ll| shows that peakons of initial width greater than a break up 
into peakons of width a under the evolution of the peakon equation 
in a periodic domain at fixed values of a = 1 and b = 2, 3. The 
emitted peakons are stable, propagate as solitary traveling waves, and 
interact elastically. Conversely, a peakon or other initial condition that 
is narrower than a will decompose into two oppositely moving trains 
of peakons and antipeakons, each of width a. 
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4.3 Pulson-Pulson interactions for b > and symmet- 
ric g 



For N = 2, the pulson dynamics in equations ([58]) and (|59| ) for b > reduces 
to 

^ = (1-^, ^ = (1-6)^, (61) 
dt y ' dqi dt y ' dq 2 



(62) 



dqi dG dq 2 dG 
dt dpi ' dt dp2 ' 

and the generating function from fl60f) is given by 

G = \{p\ +pD +pmg{qi - 92) • (63) 

(This is the Hamiltonian and the equations are canonical only for 6 = 2, 
which includes the Camassa-Holm case for which g(qi — q 2 ) = e - ' 91-92 ' gives 
the peakon solutions.) 

Conservation laws and reduction to quadrature 

Besides the total momentum 

P = Pi+P2, (64) 

the two-pulson system for b > and symmetric g also conserves a second 
quantity that is quadratic in p\ and p 2 , namely 

H = Pip 2 (l - g(qi - <?2)) 6_1 • (65) 

For a Hamiltonian system with two degrees of freedom this second con- 
servation law would be enough to ensure integrability, by Liouville's theorem. 
Even in the present case of b > without a Hamiltonian structure, this will 
be sufficient to reduce the 2-pulson system to quadratures^] 

Following the analysis for the case b = 2 and arbitrary g in Fringer and 
Holm [0, we introduce sum and difference variables as 

P = Pi + P2 , Q = qi + <?2 , v = Pi - P2 , q = qi - <? 2 • (66) 



2 When 6=1, the momenta pi and P2 are separately conserved and the problem 
immediately reduces to quadratures in q = qi — q-2 and Q = qi + q-2- 
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In these variables, the generating function (1531) becomes 



G = \p 2 (l + g(q)) + \ P 2 (l-g(q)), (67) 



and the second constant of motion (1551) becomes 



H = \{P 2 - P 2 ){l- 9 { q )) b -\ 

Likewise, the 2-pulson equations of motion transform to sum and difference 
variables as 

f =2(1-6)|| = 0, 

f = 2§ = P(1 + <?(?)), 

§ = 2(l-&)ff = W-b){P 2 -p 2 )g\q), 

I = 2 f =P(1-P(?))- 

Eliminating p 2 between the formula for H and the equation of motion for 
q yields 

{^j = P\l-g{q)f-AH(l-g{q)f~ b . (69) 

We rearrange this into the following quadrature, 

dt = -TTTTf ' Z = P\l- g( q )) 2 - AH(1 - g(q)f- b . (70) 
g'{q)VZ 



This simplifies to the quadratic Z = P 2 (l — g(qj) 2 — AH when 6 = 3. For the 
peakon case, we have g(q) = e q so that g'(q) = g(q) and the quadrature ( |70| ) 
simplifies to an elementary integral for 6 = —1,0,1,2,3. Having obtained 
q(t) from the quadrature, the momentum difference p(t) is found from fl58|) 
via the algebraic expression 



AH 
(1 "7(5) C 



P 2 = p2 -7^ 7TT7TT ' ( 71 ) 



in terms of g and the constants of motion P and iJ. Finally, the sum Q(t) is 
found by a further quadrature. The remainder of the solution for arbitrary 
6 and g closely follows Fringer and Holm for the case 6 = 2. 
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Upon writing the quantities H, P and G as 



H = cic 2 , P = ci + C2, G = -c? + -c 2 , = -P 2 



(72) 



in terms of the asymptotic speeds of the pulsons, c\ and c 2 , we find the 
relative momentum relation, 



This equation has several implications for the qualitative properties of the 
2-pulson collisions. 

Definition 4.2 Overtaking, or rear-end, pulson collisions satisfy c\c 2 > 0, 
while head-on pulson collisions satisfy c\c 2 < 0. 

The pulson order q\ < q 2 is preserved in an overtaking, or rear-end, 
collision when b > 1. This follows, as 

Proposition 4.3 (Preservation of pulson order) For overtaking, or rear- 
end, collisions when b > 1, the 2-pulson dynamics preserves the sign condi- 
tion q = qi — q 2 < 0. 

Proof. Suppose the peaks were to overlap in a collision for b > 1, thereby 
producing q = during a collision. The condition g(0) = 1 implies the 
second term in (|73|) diverges for b > 1 when the overlap occurs. However, 
this divergence would contradict p 2 > 0. □ 

Consequently, seen as a collision between two initially well-separated 
"particles" with initial speeds C\ and c 2 , the separation q(t) reaches a nonzero 
distance of closest approach q min in an overtaking, or rear-end, collision that 
may be expressed in terms of the pulse shape as, 

Corollary 4.4 (Minimum separation distance) The minimum separa- 
tion distance reachable in two-pulson collisions with c\c 2 > is given by, 




(73) 




4cic 2 



(74) 
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Proof. Set p 2 = in equation ([75[) . 



□ 



Remark 4.5 VKe s/iaZZ use result ( f7^[ ) /ater m checking the accuracy of our 
numerical simulations of these two-pulson interactions. 

Proposition 4.6 (Head-on collisions admit q — > 0) T/ie 2-pulson dynam- 
ics allows the overlap q — > w/ien b > 1 in head-on collisions. 

Proof. Because p 2 > 0, the overlap g — > implying (7 — > 1 is only possible 
in equation (|73|) with 6 > 1 for c\Ci < 0. That is, for the head-on collisions. 

□ 

Remarks about head-on collisions. For b > 1, equation (|73[) implies 
that p 2 — > 00 diverges when g — > in head-on collisions. The case b = 1 is 
regular and equation ( |73|) reduces to the constant relation p 2 = (ci — C2) 2 . For 
< 6 < 1, the quantity p 2 no longer diverges when g — > and the solution 
for the relative momentum in head-on collisions is again regular. 

4.4 Pulson-antiPulson interactions for b > 1 and sym- 
metric g 

Head-on Pulson-antiPulson collision. We consider the special case of 
completely antisymmetric pulson-antipulson collisions, for which p\ = —p2 = 
p/2 and q\ = —qi — q/2 (so that P = and Q — 0). In this case, the 
quadrature formula ([ft]) reduces to[] 

1 da ' 

^-^^Lm^' (75) 

and the second constant of motion in ( |5BD satisfies 

-AH = p 2 (l-g(q)) b - 1 . (76) 

After the collision, the pulson and antipulson separate and travel oppositely 
apart; so that asymptotically in time g(q) — > 0, p — ► 2c, and H — > — c 2 , 

3 For 6 = 3, the quadrature formula (|75| ) for the separation distance in the pulson- 
antipulson collision reduces to straight line motion, q(t) — q(to) = ±2c(t — to). 
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where c (or — c) is the asymptotic speed (and amplitude) of the pulson (or 
antipulson). Setting H = — c 2 in equation (|76|) gives a relation for the pulson- 
antipulson (p, q) phase trajectories for any kernel, 

P=± (i-Zr^- (77) 

Notice that p diverges for b > 1 (and switches branches of the square root) 
when q — > + , because g(0) = 1. In contrast, p remains constant for 6 = 1 
and it vanishes for b < 1 (and again switches branches of the square root) 
when q — > + . Note that our convention for switching branches of the square 
root allows us to keep q > throughout, so the particles retain their order. 

Remark about preservation of particle identity in collisions. The 

relative separation distance q(t) in pulson- antipulson collisions is determined 
by following a phase point along a level surface of the second constant of mo- 
tion H in the phase space with coordinates (q,p). Because H is quadratic, 
the relative momentum p has two branches on such a level surface, as in- 
dicated by the ± sign in equation (|77j) . At the pulson-antipulson collision 
point, both q — > + and either l/p — > + or p —>■ + , so following a phase 
point through a collision requires that one choose a convention for which 
branch of the level surface is taken after the collision. Taking the convention 
that p changes sign (corresponding to a "bounce"), but q does not change 
sign (so the "particles" keep their identity) is convenient, because it allows 
the phase points to be followed more easily through multiple collisions. This 
choice is also consistent with the pulson-pulson and antipulson-antipulson 
collisions. In these other "rear end" collisions, as implied by equation ([73]), 
the separation distance always remains positive and again the particles retain 
their identity. 

Theorem 4.7 (Pulson-antiPulson exact solution) The exact analytical 
solution for the pulson-antipulson collision for any b and any symmetric g 
may be written as a function of position x and the separation between the 
pulses q for any pulse shape or kernel g(x) as 

ex i 

u(x, q) = v(6-i)/2 9( x + ^/ 2 ) ~ 9(x - 9/2) , (78) 

(! -#(<?)) L J 

where c is the pulson speed at sufficiently large separation and the dynamics 



of the separation q(t) is given by the quadrature ( 75 ) with ^—AH = 2c 



35 



u(x)„ b^1 u(x), b=2 u(x), b=3 




Figure 12: Peakon-antipeakon collisions for 6 = 1,2, 3. Inviscid 6- 
family, 6 = 1, 2, 3, a = 1, initial width w — 1. 



Proof. The solution ( (57|) for the velocity u(x,t) in the head-on pulson- 
antipulson collision may be expressed in this notation as 

u(x, t) = V -g{x + q/2) - V -g(x - q/2) . (79) 

In using equation ( j77p to eliminate p this solution becomes equation ([781). ^ 



Figure [T^ shows the exact solutions for the Peakon-antiPeakon collision in 
the cases 6 = 1,6 = 2, and 6 = 3. The positive and negative peaks approach 
each other until the solution develops a negative vertical slope in finite time. 
As the separation q — > 0, the positive and negative peaks "bounce," thereby 
reversing polarity, after which they separate in opposite directions. 

4.5 Specializing Pulsons to Peakons for 6=2 and 6 = 3 

We now restrict to g(x) = e _ ' x ', the Green's function for the ID Helmholtz 
operator satisfying 

(i-spX 1 - 1 = *<*)■ w 

In this case, m = u — a 2 u xx , the pulson traveling wave solution is given by 
u(x,t) = cg(x — ct) = ce~\ x ~ ct \ has a discontinuity in derivative at its peak, 
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and is called the peakon. For 6 = 2 and 6 = 3 in the peakon case the main 
results are, 

• For 6 = 2 and g(x) = e~' x ', equation (P becomes the zero-dispersion 
limit of the integrable Camassa-Holm equation for shallow water waves 
discovered in Camassa and Holm ||. Upon restoring its linear disper- 
sion, this equation was recently proved to be a higher-order accurate 
asymptotic description of shallow water waves in Dullin et al. @ . 

• For 6 = 3 and g(x) = e - ^', equation ([!]) becomes the integrable partial 
differential equation studied in Degasperis, Holm and Hone [Q. 

• When g(x) = e"' 2- '' the N— peakon dynamics for both 6 = 2 and 6 = 3 
turns out to be integrable - as is the dynamics of the original PDE ([!]). 
The solutions of the initial value problem for (|I|) for both 6 = 2 and 
6 = 3 may be found analytically by using the Isospectral Scattering 
Transform (1ST) method. 

• The two cases 6 = 2 and 6 = 3 have quite different isospectral eigenvalue 
problems. These are discussed in Camassa and Holm [[| and in Dullin 
et al. || for the case 6 = 2, and in Degasperis, Holm and Hone 0] 
for the case 6 = 3. See also Beals, Sattinger and Smigialski JTJ] for 
a discussion of solving the inverse isospectral problem using classical 
methods for the case 6 = 2. 



5 Peakons of width a for arbitrary b 

When g = e~^ x ^ a , we may invert the velocity-momentum relation u = g * m 
by using the Green's function expression (^) with the Helmholtz operator 
to find m = u — a 2 u xx . Hence, equation (fj) may be rearranged into the local 
momentum conservation law, 

d ( . & - 1 2 b- 1 2 2\ / a n 
mt = ~ dx \ ~^T U 2~ Q U V ' ^ ' 

This conservation law for peakons may also be rewritten in convection 
form: 

Ut + uu x = — t x with (1 — a o x )r = -u H — a u x . (82) 
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The two forms (|5T| ) and (^) of the b-family of equations ([!]) suggest that 
values b = 0,1,3 are special. These values of b are natural candidates for 
boundaries, or bifurcation points for changes in solution behavior. 

Equation (|S2|) describes peakons of shape g(x) = e~^ x ^ a . This peakon 
equation will form the basis of the rest of our study. 

5.1 Slope dynamics for Peakons: inflection points and 
the steepening lemma when 1 < b < 3 

We shall consider solution dynamics of equation fl82|) in the peakon case sat- 
isfying (||^), or equivalently, equation ([I]) with g(x) = e~^ x ^ a , which satisfies 

(1 - a 2 d 2 x )e-W /a = 2a5(x) . (83) 

For this case, and with vanishing boundary conditions at spatial infinity, 
equations (1521) and (1531) imply the peakon equation on the real line, 



oc 



/// ////, = — / e ~\ x ~y \l a (buuy + (3 - b)a 2 u y u yy ) dy. (84) 



2a 

Taking the x— derivative gives the equation for the slope u x (x,t) 

1 d ro ° 



U xt + UU XX + III = 

2a dx 



I e~ lx ~ y l/a (buu y + (3 - b)a\u yy ) dy 

J —oo 

^ [ e~ lx ~ y l/a sgn (x - y) (buu y + (3 - b)a 2 u y u yy \ dy . (85) 

^ J — oo 



2a 

We shall use these expressions to prove the following. 

Proposition 5.1 (Peakon Steepening Lemma) For b in the range 1 < 
b < 3 a sufficiently negative slope at an inflection point of u will become 



vertical in finite time under the dynamics of the peakon equation (84)- 



Proof. Following Camassa and Holm 0, we shall consider the evolution 
of the slope u x at an inflection point x = x(t). Define the slope at the 
inflection point as s(t) = u x (x(t),t) and note that u xx (x(t),t) = 0. Then 
equation (|85|) yields the following evolution equation for s(t) 

^ t + s2 = ^2 /e-'^^sgn (*(*) - y)\§- y (bu 2 + (3 - b)a 2 u 2 y ) dy (86) 
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Integrating by parts using the definition \y\ = ysgn(y), so that d\y\/dy 
sgn (y) + 2y8(y), and recalling that yS(y) = 0, gives 



Hence, in the range < 6 < 3 the last term is negative and we have the slope 
inequality, 

We suppose the solution satisfies (bu 2 /a 2 ) < M for some constant M.f\ Then, 

ds (b - 1\ , M „ 

— < — ( ) s 2 + — for < b < 3 . 

dt ~ V 2 / 2 ~- 

Consequently, if 6 > 1, 



dcoth-^X) < y/M for X = \ h —^s. (89) 



1 _X 2 - ' "" ~" V M 

This implies, for s < —\/M initially negative, that 



s < VMcoth(<7 + y-jj^Y*) for 1<^<3, (90) 

where the dimensionless integration constant a < determines the initial 
slope, which is negative. Under these circumstances, the slope at the inflec- 
tion point must become vertical by time t = —2a/ a/ M(b — 1). 

□ 



4 If this inequality is violated, we have another type of singularity. However, for 6 = 2, 
the constant M can be estimated by using a Sobolev inequality In fact, M = AH\(u)/a 2 
because for this case we have 

1 f°° 

max[w 2 (x, t)] < — / (u 2 + a 2 ul) dx = 2H 1 = const for 6 = 2. 
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Remarks for 1 < b < 3. 

• If the initial condition is antisymmetric, for 1 < b < 3, then the inflec- 
tion point at it = is fixed and dx/dt = 0, due to the mirror reflection 
symmetry (u,x) — > (—u, —x) admitted by equation fl34|). In this case, 
M = and equation (]88[) implies 

4.1* - »w<^(^), OK 

so verticality s = — oo will develop in finite time, regardless of how 
small the initial slope |s(0)|, provided it is negative, s(0) < 0, as in 



figure 12. If the initial slope is positive, then under this evolution it 



will relax to zero from above. 

Consequently, traveling wave solutions of ( |84T ) cannot have the usual 
sech-like shape for solitons because inflection points with sufficiently 
negative slope can produce unsteady changes in the shape of the solu- 
tion profile. 



5.2 Cases < b < 1 

In the range < b < 1, we have from (^) that 

/oo roo 
|m| 1//f) <ix = / |m | 1//f) dx , where mo(x) = m(x, 0) . (92) 
-oo J — oo 

This conservation law implies an elliptic regularity estimate showing that the 
slope s = u x is always bounded under the dynamics of the peakon equation 
(|84p. See Holm and Titi [2002] for a proof of this result and more discussion 
of its implications. 



6 Adding viscosity to peakon dynamics 

In the remainder of this paper, we shall restrict our attention to the peakon 
case g(x) = e~\ x \l a with length scale a, and investigate the fate of the peakon 
solutions when viscosity is introduced for given values of b and a. For pur- 
poses of comparison with previous results in the literature, we shall also 
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extend equation (]!]) to a new family of equations that includes the Burg- 
ers equation by introducing two additional real parameters. These are the 
viscosity v and a multiplier j3 for the stress, or pressure gradient. 

First, we shall introduce constant viscosity v > into (P to form the 
viscous b-family of equations for the peakon case g(x) = e~^ x ^ a , as follows, 

m t + um x + bu x m = vm xx , with m = u — a 2 u xx . (93) 
convection stretching viscosity 

As in equation (^), this equation with viscosity may be expressed solely in 
terms of the velocity u(x,t) as 

u t + (b + l)uu x - vu xx = a 2 (u xxt + uu xxx + bu x u xx - vu xxxx ) (94) 

= a 2 d x (u xt + uu xx - vu xxx + ^—u 2 x 

o„ 2 / b — 3 2 N 
= a d x \u t + uu x - vu xx H — u% 

Thus, the nonlinear steepening term increases with b as (b + l)uu x . When 
a — » the previous equation reduces to 

u t + (b + l)uu x - vu xx = , (95) 

and one then recovers the usual Burgers equation either by rescaling di- 
mensions, or by setting b = 0. For b = 2, equation (|93| ) is the one-dimensional 
version of the three-dimensional Navier- Stokes- alpha model for turbulence 

The viscous b-family of peakon equations (^) may be rearranged into two 
other equivalent forms that are convenient for making the second extension 
of a stress multiplier. These are either its equivalent conservative form, 

d ( b — 1 2 b — 1 22 \ .. 
m t = - — [mu H —u — a u x J + vm xx , (96) 

or its equivalent convective form, 

(1 - a 2 d 2 x ) (u t + uu x - vu xx ) = - d x (^u 2 + ^-^a 2 u^ . (97) 
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Stress multiplier (3. Next, we shall introduce a stress multiplier (3 as a 
second parameter that for (3^1 deforms the convective form of the viscous 
b-family of equations fl9~T|) into the following family of Burgers-like equations 
with four parameters b, a, v and (3, 

b 3 b 

u t + uu x - vu xx = - (3 t x with (1 — a 2 d 2 )T = -u 2 H — a 2 u 2 x . (98) 

When (3 = 0, the Burgers— a(3 equation (p8|) recovers the usual Burgers 
equation. When (3 = 1, equation (|98|) recovers the viscous b-family of peakon 
equations (p3|). 

We shall seek solutions of the Burgers— a/3 equation (|98|) , either on the 
real line and vanishing at spatial infinity, or in a periodic domain, for various 
values of its four parameters b, a, v and [3. Under these boundary conditions, 
when /3 — *■ 1, equation ( |9"8| ) recovers the convective form (£J7p of the viscous 
b-family for peakons with g(x) = e~' x ' / ' a . Thus, the viscous b-family of 
equations (|93| - p7D deforms into the Burgers— a(3 equation (|98|) when (3^1 
and the Burgers— a/3 equation (|98| ) reduces to the usual Burgers equation 
when (3 = 0. We shall be interested in the effects of the four parameters b, 
a, v and (3 on the solutions of the Burgers— a/3 equation d98|). We shall be 
interested especially in the fate of the peakon solutions upon introducing the 
parameters v and (3 so as to retain H\ control of the velocity. As we shall 
see, such control requires a special relation between the parameters b and j3, 
namely, (3 — b)j3 = 1. 

6.1 Burgers— af3 equation: analytical estimates 

Proposition 6.1 (H^ control of the velocity) The Burgers— a(3 equation 
( p^ ) controls the a— weighted H 1 norm of the velocity, 

/oo 
(u 2 + a 2 u 2 x ) dx 
■oo 

for a 2 7^ 0, provided (3 — b)(3 = 1. 
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Proof. The spatial derivative of the Burgers— a/3 equation ([53) yields the 
dynamics for the slope s 



U t + UU X - VU XX = —@T X , 

s t + us x + s 2 — vs xx = — (3 t xx , with s = u x 



9 b 9 3 b o 9 

-a = —u H a u„ — t . 



In turn, these slope dynamics equations imply the following evolution of the 
a— weighted H 1 density, cf. equation ([86|), 

d (\ 9 a 2 9 \ d f 1 / b/3\ , _ a 2 

4 + ^((3-6)/3-l) S 3 . 



Thus, provided 

(3 - b)P = 1 , 

the last term will vanish. Under this condition, for periodic or vanishing 
boundary conditions the a— weighted H 1 norm 

/oo 
(u 2 + a 2 w 2 ) dx 
-oo 

will decay monotonically under the Burgers— a/3 dynamics for a 2 ^ 0. 

□ 



Remarks. 

• When v — > in the Burgers— a/3 equation, the a— weighted if 1 norm 
is conserved for (3 — 6)/3 = 1. This relation cannot be satisfied for 
6 = 3. Thus, the proof of decay of the a— weighted H 1 norm under the 
Burgers— a/3 dynamics is inconclusive for v ^ when 6 = 3. However, 
one can expect on physical grounds that this norm will also decay for 
b = 3 if v is sufficiently large. 

• We shall restrict our remaining considerations to those values of b and 
P for which the a— weighted H 1 norm is bounded, or decays mono- 
tonically. In one dimension, this control of the a— weighted H 1 norm 
implies the solution for the velocity will be continuous. 
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• Namely, we shall consider the following cases with (3 — b)(3 = 1 
(b=0, (3 = 1/3), (b=l, (3 = 1/2) and (b=2, (3 = 1). 

Proposition 6.2 (Burgers— af3 Steepening Lemma) Forb and (3 in the 
range (3—b)j3 < 2 a sufficiently negative slope at an inflection point of velocity 
u will become vertical in finite time under the dynamics of the Burgers— a/3 



equation ( \98j) with v = 0. 



Proof. The proof follows that for the Peakon Steepening Lemma |5J] and 
uses the slope equation following from Burgers— a/3 equation fl98|) with v — 
that corresponds to fl85|) for the Peakons, modified to include f3, 



u xt + uu xx + u 2 x = - T^j^; J e lx yl/a (buu y + (3 - b)a 2 u y u yy ^ dy 



f3 d 



2a 



f e -\x- y |/« ggn ^ _ y ^ f buUy + (3 _ b)a 2 u y u yy \ dy . (99) 

J —oo 

Equation (|99| ) yields the inviscid Burgers— a (3 evolution of the slope s(t) = 
u x (x(t),t) at an inflection point x = x(t) as 

|<-(2-(3-6)/?)^ + ^ for 0<6<3, (100) 

This holds provided we assume the solution satisfies (bu 2 /a 2 ) < M for some 
constant M. Consequently, if 2 — (3 — b) (3 > 0, we have 



dX 



dcoth-^X) < VM for X = J- — t^-t^-s. (101) 



i - x 2 v / - y Mp 

For s < —y/M initially negative and j3 > 0, this implies, 



s < v / Mcoth(a+ J- — ^ 1) for < 3 - | < b < 3 , (102) 

where the dimensionless integration constant a < determines the initial 
slope, which is negative. Under these circumstances, provided the inflection 
point continues to exist, its negative slope must become vertical by time 
t - =2£. I M P 

L M V 2-(3-fc)/T 

□ 
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Corollary 6.3 (Inviscid Burgers— a/5 shocks) Solutions of the inviscid 
Burgers— aj3 equation with v = that remain continuous in velocity 
must develop negative vertical slope in finite time. 



Proof. According to Proposition |6.1|, continuity of the velocity and, hence, 



control of the H l norm ||w||#i requires that (3 — b)f3 = 1. This is in the 
parameter range where Proposition |6.2| applies. Consequently, verticality 
will form at an inflection point of negative slope under the dynamics of the 
inviscid Burgers— a/3 equation (BSf) with v = for (3 — b)j3 < 2. □ 



Remark 6.4 Hence, to remain continuous without viscosity, the solution 
of the inviscid Burgers— a(3 equation must either develop verticality at an 
inflection point of negative slope, or it must evolve to eliminate such points 
entirely. 

6.2 Burgers— a/3 traveling waves for /?(3 — b) = 1 &; v = 



For v = 0, the Burgers— a/3 equation (|98|) has traveling waves given by 



(u-c)u' + Pt' = and r - aV = h -u 2 + ^-^ a 2 ( u ') 2 (103) 
which yields after one integration 

cu + pr = K, (104) 



u 2 



2 

where K is the first integral. Consequently, we find 

,2 



1 / U \ 

r-<rr- = -[K + cu- — + a 2 ({u-c)u"+{u') 2 )) . (105) 



. 2_ // _ 



The second equation in ( |103| ) integrates for the special case of (5 (3 — b) = 1, 
2Ku + cu 2 - pu s + a 2 (u - c)(u ') 2 = 2H . (106) 
For the special case K = = H this becomes 

a 2 {u-c){u') 2 = ({3u-c)u 2 for /3(3 - b) = 1 , (107) 
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and we recover the peakon solution u(z) = ce - ' 2 '/" for (3 = 1. In the general 
case that K ^ and H ^ 0, we rearrange equation ( |106|) into the following 
quadrature for inviscid Burgers— a/3 traveling waves, 

dz (u — c) 1 ' 2 du „, „ . . 

± — = i - for /3(3 - b) = 1 & v = . (108) 



a 



2H - 2Ku - cu 2 + pu 3 



In what follows, we shall consider the cases (6 = 0,/? = 1/3), (b = 1, (3 = 1/2) 
and (6 = 2, /? = 1) when v ^ 0. 

7 The fate of the peakons under (1) adding 
viscosity and (2) Burgers— a(3 evolution 

7.1 The fate of peakons under adding viscosity 

The following set of four figures shows the effects on the initial value problem 



for the viscous b-equation (|93|) of varying a and b at fixed viscosity for an 
initial velocity distribution given by a peakon of width w = 5 and initial 
height U ~ 0.1. The parameter b takes the values b = 0, 1, 2, 3. In these four 
figures, the resolution is 2 13 points on a domain size of 200 with viscosity 
v = 0.005. This corresponds to a grid-scale Reynolds number of Re& x = 
UAx/i/ = 0(1) for velocity U — 0.1. The pair of figures after these four then 
shows the effects on the same problem of increasing viscosity v at fixed a for 
6 = 2 and 6 = 3. 



Figure [T3| shows three plots of the evolution of the velocity profile under the 
viscous b-equation (|93|) of an initial peakon of width five, as a function 
of increasing a = 1/4, 1,4 at fixed viscosity v = 0.005 for 6 = 0. The 
peakon leans to the right and develops a Burgers-like triangular shock, 
or ramp and cliff, whose width increases and peak height decreases as 
a increases. These three plots show no discernable differences for 6 = 
as the viscosity is decreased to v = 10 -6 . Hence the width of the cliff 
in the ramp and cliff structure for 6 = is set by the value of a in this 
range of parameters. 



Figure [14] shows three plots of the same type of evolution from a peakon 



initial condition of width w = 5, as a is varied for 6 = 1. The front of 
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the ramp and cliff structure propagates faster and is sharper for 6=1 
than for 6 = when a = 1/4 and a = 1. This increase in speed 
appears to occur because the coefficient increases in the steepening 
term (6 + l)uu x in equation (P4"|). A nascent peakon begins to form 
close behind the front at the top of the ramp, then eventually gets 
absorbed into the ramp and cliff. For a = 4, however, this nascent 
peakon forms more completely and nearly escapes. 



Figure ITS] again shows three plots of the evolution from a peakon initial 



condition of width w = 5, as a is varied, this time for 6 = 2. The ramp 
and cliff structure is faster for 6 = 2 than for 6 = 1 when a = 1/4. 
When a = 1 a series of three nascent peakons forms close behind the 
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U(X) b=2 o=1/4,1,4 T=0,500,1000. 1500,2000,2500 

0.1 5i 1 i 1 p r 1 1 1 t 




X 

Figure 15: Effect of increasing a for 6 = 2. Viscous 6-family, 6 = 2, 
a = 1/4, 1,4, v = 0.005, initial width w = 5. 



front, then overtakes the ramp and cliff structure and slightly affects its 
propagation before eventually being absorbed. For a = 4, however, the 
initial peakon simply propagates and decays under viscosity, although 
it is slightly rounded at the top. 



Figure |16| also shows three plots of the evolution from a peakon initial 
condition of width w = 5, as a is varied, this time for 6 = 3. The 
ramp and cliff structure moves faster yet, and a single nascent peakon 
appears just behind the front already for a = 1/4. When a = 1, a series 
of three nascent peakons forms initially close behind the front and they 
nearly escape before being slowed by viscosity. The leading peakon 
decays and slows due to viscosity. Then the following ones overtake 
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U(X) b=3 o=1/4,1,4 T=0,500,1000. 1500,2000,2500 

0.1 5i 1 i 1 p r 1 1 1 t 




X 



Figure 16: Effect of increasing a for 6 = 3. Viscous 6-family, 6 = 3, 
a = 1/4, 1,4, v = 0.005, initial width w = 5. 



and collide with the ones ahead as the ramp and cliff structure forms. 
These collisions occur at higher relative velocity for 6 = 3 than for 6 = 2 
and they significantly affect the propagation and eventual formation of 
the ramp and cliff. In contrast, for a = 4, the initial peakon keeps its 
integrity and simply propagates rightward and decays under viscosity. 
The propagating peakon for a = 4 at this viscosity decays more slowly 
and is much sharper at the top for 6 = 3 than for 6 = 2. 



Remark 7.1 (Exchange of stability) To see the exchange of stability be- 
tween the ramp /cliff structure and the peakon as b changes, we perform the 
following numerical experiment. First, we run the viscous b- equation ( ffi^ j 
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with 6 = 0, a = 1, v = 10 , and an initial peakon of width w = 5. As 
we see in Figures [7?| and |7^, this evolves into the ramp and cliff formation 
even for nearly zero viscosity. Once the final ramp/cliff state is formed, we 
then use it as the new initial condition for equation ffijj ) with either b = 2 or 
b = 3. The new evolution breaks the ramp/cliff structure into peakons and 
the new final state is a rightward moving train of peakons ordered by height. 

For Figures and \2Q, we ran the same numerical experiment, this time 
with a value a = 5 equal to the width of the initial peakon. The initial peakon 
"borrows from the negative" to form a ramp, which is not quite antisymmetric 
because the total area of the initial peakon must be preserved. At time T = 150 
we switch to b = 2 (top plot) or b = 3 (bottom plot), and again observe a 
train of stable peakons emerging from the now-unstable ramp. 

Finally, for Figures and \2^, we again run the numerical experiment 
with a = 1 and an initial peakon width w = 5, but this time changing to 
b = —2 or b = — 3 after the ramp has formed. The new evolution breaks the 
ramp/cliff structure into leftons like those in Figures^ and^. 



Remark 7.2 (Increasing viscosity) The effect of increasing viscosity on 
the evolution of the peakon initial condition can be estimated from the a— scale 
Reynolds number defined by, 

Re a = Uajv = (a/ Ax)ReAx ■ 

For (a/ Ax) = 40, U = 0.1 and increasing viscosity v, the Reynolds numbers 
Re^x and Re a decrease as 

Re Ax ^ 2, 0.2, 0.02 and Re a =~ 80, 8, 0.8 for v = 0.01, 0.1, 1.0 . 

Perhaps not surprisingly, when Re a = 0(1) the viscosity will diffuse through 
the initial peakon before it can fully form. Figures [|3] and |^| show that this 
effect increases as Re a decreases. 



7.2 The fate of peakons under Burgers— af3 evolution 



Figures £25] and show the effects on the peakon initial value problem for 
the Burgers— a/3 evolution of varying a and b with (3 — b)f3 = 1 at constant 
viscosity. We shall consider the following cases with (3 — b) (3 = 1: 
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b=0 ta T=2500, then b=2 to T=5000 




Figure 17: Exchange of stability between ramps and peakons for 
b = 0,2,3, when width > a. Viscous 6-family, b = — > 2,3, a = 1, 
v = 10~ 5 , initial width w — 5. 
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T=0, T=2500 (switch from b=0 to b=2). T= 500.0 




20 40 60 



80 100 120 .140 100 ISO 
X 



200 



T=0, T-2500 (switch from b=0 to b=3) a T-5000 




l 40 160 



Figure 18: Exchange of stability between ramps and peakons for 
b = 0, 2, 3, when width > a: profiles. Viscous 6-family, 6 = 0^2,3, 
a — 1, v — 10~ 5 , initial width w — 5. 
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b=0 ta T=1 50. then b=2 to T=400 




-200 -150 -100 -50 50 100 150 200 

X 



Figure 19: Exchange of stability between ramps and peakons for 
6 = 0,2, 3, when width = a. Viscous 6-family, b = — > 2, 3, a = 5, 
v = 10~ 5 , initial width w — 5. 
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T=0, T=1 50 (switch from b=0 to b=2). T=300 




T=0, td 50 (switch from b=0 to b=3) s T=300 
3i 1 1 1 1 i i 1 1 



3 




Figure 20: Exchange of stability between ramps and peakons for 
b = 0, 2, 3, when width = a: profiles. Viscous 6-family, 6 = 0^2,3, 
a — 5, v — 10~ 5 , initial width w — 5. 
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b=0 to T=2500, then b=-2 to T=7500 




Figure 21: Exchange of stability between ramps and leftons for 
6 = 0, —2, —3, when width > a. Viscous Warmly, b = ^ —2, —3, a — 1, 
v = 10~ 5 , initial width w — 5. 
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T=0, T=2500 (switch from b=0 to b=-2), T=7500 




T=Q S T=2500 (switch from b=0 to fe=-3), T=7500 




Figure 22: Exchange of stability between ramps and leftons for b = 
0, —2, —3, when width > a: profiles. Viscous 6-family, b = — > —2, —3, 

a — 1, v — 10~ 5 , initial width w — 5. 
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u(x) b=3 v=1/1 00, 1/10,1 T=0.500 ; 1 000,1 500 .2000,2500 
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ure 24: Effect of increasing viscosity for 6 = 3. Viscous 6-family, 
3, a = 1, v = 1/100, 1/10, 1, initial width w = 5. 
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b = 0, = 1/3, v = 0.005, a = 1/4, 1, 4, and 
b = 1, p = 1/2, v = 0.005, a = 1/4, 1, 4. 



Remark 7.3 (Lowering P has little effect on the ramp/cliff) Lowering 
P to follow (3 — b)P — 1 instead of keeping P = 1 /ias Zztt/e effect on i/ie (ie- 
velopment of the ramp/cliff solution for 6 = and 6=1. Lowering P for 
these cases only makes the activity slightly less lively at the front for (6 = 0, 
P = 1/3) and (b = 1, P = 1/2) than for the corresponding cases ofb = and 
6=1 with P = 1 in Figures [JJ| and [T^J. This lessened activity at lower P can 
only be discerned in the solution for the largest value a = 4. The remaining 
case (b=2, P = 1) recovers the viscous b-equation (jP^) for 6 = 2 in Figure 
1731 , in which the larger b produces much livelier steepening and, hence, more 
activity at the front of the rightward moving pulses. 



8 Numerical results for peakon scattering and 
initial value problems 

We shall begin by summarizing the results in the figures given earlier, and 
then describe the numerical methods used in producing them and discuss 
some of the ways we verified and validated the results. 

8.1 Peakon initial value problems 

8.1.1 Inviscid b-family of equations 

Ramps and cliffs for 6 = 0. Figure [I] shows the formation of a ramp and 
cliff pattern for 6 = 0, a = 1, and a set of Gaussian initial conditions 
of increasing width w = 2.5, 5, 10. 

Peakons for 6 = 2, 3. Figures |2| and |3| show the formation of peakons for 
6 = 2 and 6 = 3, for a = 1 and a set of Gaussian initial conditions of 
increasing width w = 2.5, 5, 10. 

Ramps and cliffs for 6 = —1/2. Figure f| shows the formation of a ramp 
and cliff pattern for 6 = —1/2, a = 1, and a set of Gaussian initial 
conditions of increasing width w = 10, 15, 20. 
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Figure 25: Effect of increasing a when (3 — b)f3 = 1, for 6 = and 
(3 = 1/3. Burgers-a/3, b = 0, a = 1/4, 1, 4, (5 = 1/3, v = 0.005, initial width 
w — 5. 
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Stationary solutions for b < — 1. Figure [| shows an essentially station- 
ary solution with a slight rightward drift and leaning slightly to the 
right due to nonlinear curvature terms with higher order derivatives 
in equation fl4~9|), for a = 1 and a set of Gaussian initial conditions of 
increasing width w = 10, 15, 20. For the same a and same set of initial 
conditions, Figures ^| and [7j show the emergence of leftons. 

Figure § shows the leftons at time T = 2500 for the b = —2 case, 
versus the analytical u(x) ~ sech 2 (a;/(2a)) from equation ((52|), and for 
the b = —3 case, versus the analytical u{x) ~ sech(x/a). 

Peakons of width a for b = 2, 3. Figures [5] and [lTj show 2-peakon and 
3-peakon interactions for 6 = 2 and for 6 = 3, beginning with initial 
peakons of width w = a = 5. 



Peakons of width > a for 6 = 2,3. Figure 11 shows the emergence of 
peakons of width a = 1 when we begin with peakons of width w = 5 
greater than a, for 6 = 2 and 6 = 3. 

Peakon-antipeakon collisions for 6 = 1,2, 3. Figure [1^ shows the dy- 
namics of a peakon-antipeakon collision for 6=1,6 = 2, and 6 = 3, for 
a = 1, at four successive times. 

8.1.2 Viscous b-family of equations 
Effect of a for 6 = 0, 1,2,3. Figures O - [TBI show the evolution of an 



initial peakon of width w = 5 as a function of increasing a = 1/4, 1,4 
at fixed viscosity v = 0.005, for 6 = 0, 6 = 1, 6 = 2, and 6 = 3. 



Exchange of stability between ramps and peakons. Figures |lj] and [TJ 
show the exchange of stability between ramps and peakons suggested 
in the previous four figures, with a = 1 and an initial peakon of width 
w — 5, but this time with a very small viscosity v = 10~ 5 so that 
the peakons, when stable, do not noticeably decay. The exchange of 
stability occurs when we switch from 6 = to 6 = 2 or 6 = 3. Figures 
and |2pj again show the exchange of stability, this time using a = 5 
so that the initial peakon has width a. 



Exchange of stability between ramps and leftons. Figures |21] and ^2 



show the exchange of stability as in Figures ITT] and ITS, but we switch 
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to b = —2 or b = —3 instead, and see the emergence of stable leftons 
from the ramp. 



Effect of viscosity for b = 2,3. Figures El and 121 show the effect of 



increasing viscosity v = 0.01, 0.1, 1 on the evolution of an initial peakon 
of width w = 5, with a = 1 and b = 2 or b = 3. 



8.1.3 Burgers-a/3 equation 
Effect of a when (3 — 6)/? = 1. Figures 123 and 126] show the effect of 



increasing a = 1/4, 1,4 on the evolution of an initial peakon of width 
w = 5, for fixed v = 0.005 and two sets of values for b and (3 for which 
(3 - b)(3 = 1: b = 0, 13 = 1/3 for the first figure, and b = 1, (3 = 1/2 for 
the second figure. 

8.2 Description of our numerical methods 



For our numerical runs we advanced equations fl82"D, (97), and ( j9~8|) with an 
explicit, variable timestep fourth/fifth order Runge-Kutta-Fehlberg (RKF45) 
predictor/corrector. We selected the timestep for numerical stability by trial 
and error, while our code selected the timestep for numerical accuracy (not 
to exceed the timestep for numerical stability) according to the well-known 
formula from numerical analysis, 

f Ah I \ 1 / p 

hi = 7^-1 ,, ' ,, • (109) 



This is used in the following way. At step i of the calculation, we know the 
predicted solution Ui, the corrected solution u h and the previous timestep 
hi-i. The predictor's order of accuracy is p = 4, while the corrector's order 
of accuracy is p+ 1. A new timestep hi is chosen from ( |109| ) based on the old 



timestep hi-\ and the norm of the difference between the current predicted 
and corrected solutions. We used a very strict relative error tolerance per 
timestep, e = 10~ 8 , a safety factor 7 = 0.9, and an L 2 norm || ■ || 2 . 

We computed spatial derivatives using 4th order finite differences, gen- 
erally at resolutions of 2 13 = 8192 or 2 14 = 16384 zones. To invert the 
Helmholtz operator in transforming between m(x,t) and u(x,t), we con- 
volved m(x, t) with the Green's function in Fourier space. When the nu- 
merical approximation of the nonlinear terms had aliasing errors in the high 
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wavenumbers, we applied the following high pass filtered artificial viscosity, 




if < k < f , 



"(*) = <¥(*-#) if f < fc <¥> ( n °) 

if ^ < k < N, 

where 5 = 0.01 for the present simulations. iV is one-half the number of 
zones, because for each zone we have both a Fourier sine coefficient and a 
Fourier cosine coefficient. 

The quality of the numerical convergence may be checked analytically in 
the case of rear-end two-pulson collisions, for which equation ([74]) in Corollary 
fO| yields 

( \ (\ i ( 4ciC2 V /(b ~ 1} run 

g{q m in)=g{q) = 1 - ( 7 - , 2 ) • HI 

For peakons with 6 = 2 and = e - '^/", this formula gives the minimum 
separation, 

gmin = -2aln( > 0. (112) 

V Ci + c 2 / 

When ci = 1, C2 = 1/2, and a = 5, as in figure |9|, this formula implies q m i n = 
10 In 3 = 10.9861. Our numerical results with the resolution of 2 14 zones 
yield q m i n = 11.0049. The very small discrepancy, less than 0.2%, occurs 
largely because our numerical measurement of q m m is obtained by examining 
the peakon positions at each internal timestep in the code, while the code's 
time discretization effectively means we're unlikely to land exactly on the 
time at which the minimum separation occurs. The code's true accuracy is 
thus better than the above measure indicates, because the intermediate steps 
involved in advancing the solution from one discrete time to the next with 
an RKF45 method cancel the higher-order discretization errors. 

Likewise, for peakons with 6 = 3 and g(x) = e~' x '/ a , formula ( 111 ) gives 
the minimum separation, 



q min = -a\n( 1 - T ^^— ) >0. (113) 
^ (ci + c 2 )/2/ 

When Ci = 1, c 2 = 1/2, and a = 5, as in figure [l(| this formula implies 
q-min = 5 In (3/(3 — a/8)) = 14.3068. This time our numerical results yield 
q m in = 14.2924, a discrepancy of only 0.1%. 
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Of course, the two-body collision is rather simple compared to the plethora 
of other multi-wave dynamics that occurs in this problem. For this reason, we 
also checked the convergence of our numerical algorithms by verifying that 
the relative phases of the peakons in the various figures remained invari- 
ant under grid refinement. Moreover, the integrity of the waveforms in our 
figures attests to the convergence of the numerical algorithm - after scores 
of collisions, the waveforms given by the Green's function for each case are 
still extremely well preserved. The preservation of these soliton waveforms 
after so many collisions would not have occurred unless the numerics had 
converged well. 

9 Conclusions 

Equation ([]]) introduced a new family of reversible, parity invariant, evolu- 
tionary 1 + 1 PDEs describing motion by active transport 

m t + um x + bu x m = 0, with u = g*m. (H4) 
convection stretching 

We analyzed the transformation properties and conservation laws of this 
family of equations, which led us to choose g to be an even function. Then 
we classified its traveling waves, identified the bifurcations of its traveling 
wave solutions as a function of the balance parameter b and for some choices 
of the convolution kernel g(x) we studied its particle-like solutions and their 
interactions when b > 1. These were obtained by superposing N traveling 
wave solutions u(x,t) = cg(x — ct) as 

N N 

u(x, t) = y^ y pj(t)g(x - qj{t)) and m(x, t) = y^j^Sjx - , (115) 

i=i i=i 

for any real constant b and u = g * m, in which the function g is even 
g(— x) = g(x), so that g'(0) = 0, and is bounded, so we may set g(0) = 1. 

Following 0, we call these solutions "pulsons." We have shown that for 
any b > 1, once they are initialized on their invariant manifold (which may 
be finite dimensional), the pulsons undergo particle-like dynamics in terms of 
the moduli variables Pi(t) and qi(t), with i = 1, . . . , N. The pulson dynam- 
ics we studied for b > 1 in this framework on a finite-dimensional invariant 
manifold displayed all of the classical soliton interaction behavior for pulsons 
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found in JT] for the case 6 = 2. This behavior included pairwise elastic scat- 
tering of pulsons, dominance of the initial value problem by confined pulses 
and asymptotic sorting according to height - all without requiring complete 
integrability. Thus, the "emergent pattern" for b > 1 in the nonlinear evolu- 
tion governed by the active transport equation (H) was the rightward moving 
pulson train, ordered by height. Thus, the moduli variables Pi(t) and qi(t) 
are collective coordinates on an invariant manifold for the PDE motion gov- 
erned by equation ([]]). Once initialized for b > 1, these collective degrees of 
freedom persist and emerge as a train of stable pulses, arranged in order of 
their heights, that then undergo particle-like collisions. 

In contrast, the emergent pattern in the Burgers parameter region < 
b < 1 is the classic ramp/cliff structure as in Figure O. That the behavior 



should depend on the value of b is clear from the velocity form of equation 
([]]) written in 



u t + (b + l)uu x - uu xx = a 2 (u xxt + uu xxx + bu x u xx - vu xxxx ) (116) 

= a 2 d x (u xt + uu xx - vu xxx + — — ^ ul 

2„,/ b — 3 2 N 
= a d x yu t + uu x - vu xx H — u x 

Thus, nonlinear terms in this equation change sign at four integer values of 
the parameter b. Nonlinear a 2 — terms change sign when b = 0,1,3. Also, 
the nonlinear steepening term increases with b as (b + l)uu x . So this term 
changes sign when b = — 1. In the parameter regime b > — 1 (resp. b < — 1) 
the solutions of equation (|l|) move rightward (resp. leftward), provided the 
terms on the right hand side of equation ( [L16| ) are sufficiently small. 



Three regions of b. We found that the solution behavior for equation ([]]) 
changes its character near the boundaries of the following three regions in 
the balance parameter b. 

(Bl) In the stable pulson region b > 1, the Steepening Lemma for peakons 
proven for 1 < b < 3 in Proposition [5j] allows inflection points with 
negative slopes to escape verticality by producing a jump in spatial 
derivative at the peak of a traveling wave that eliminates the inflection 
points altogether. Pulson behavior dominates this region, although 
ramps of positive slope are also seen to coexist with the pulsons. When 
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b < 1 we found the solution behavior of the active transport equation 
(P changed its character and excluded the pulsons entirely. 

(B2) In the Burgers region < b < 1, the l}l h norm of the variable m 
is controlled^] and the solution behavior is dominated by ramps and 
cliffs, as for the usual Burgers equation. Similar ramp/cliff solution 
properties hold for the region — 1 < b < 0, for which the L 1 ' 6 norm of 
the variable l/|m| is controlled. At the boundary of the latter region, 
for b = — 1, the active transport equation ([[]) admits stationary plane 
waves as exact nonlinear solutions. 

(B3) In the steady pulse region b < — 1, pulse trains form that move leftward 
from a positive velocity initial condition (instead of moving rightward, 
as for b > —1). These pulse trains seem to approach a steady state. 



Effects of viscosity. Almost any numerical investigation will introduce 
some viscosity or other dissipation. Consequently, we studied the fate of the 
peakons when viscosity was added to the b-family in equation fl93|). Viscous 
solutions of equation (P3I) for the peakon case g(x) = e~\ x \l a with a = 1 



were studied in each of the three solution regions (B1)-(B3). In the Burgers 
region (B2) near b = we focused on the shock- capturing properties of the 
solutions of equation ([]]) and this family of equations was extended for (3 ^ 1 
to the Burgers— a (3 equation (|98|) , 

b 3 b 

u t + uu x — vu xx = — (3 t x with (1 — a 2 d x )r = —u 2 H — a 2 u x . (117) 

According to Proposition |6.1| , the Burgers— a/3 equation ( |1 1 7] ) controls the 
a— weighted H 1 norm of the velocity for a 2 ^ 0, provided (3 — b)f3 = 1. This 
analytical property guided our study of this new equation by identifying 
a class of equations for which a priori estimates guarantee continuity of the 
solution u(x, t). The shock-capturing properties of the Burgers— a(3 equation 
(|117|) and its a — > oo limit will be reported in a later paper fTCf, 
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